Constraints on the induced gravitational wave background 
from primordial black holes 



Edgar Bugae\0 and Peter Klimajll 

Institute for Nuclear Research, Russian Academy of Sciences, 
60th October Anniversary Prospect 7a, 117312 Moscow, Russia 

We perform a consistent calculation of primordial black hole (PBH) mass spectrum and second- 
order induced gravitational wave (GW) background produced from primordial scalar perturbations 
in radiation era of the early Universe. It is shown that the maximal amplitudes of the second or- 
der GW spectrum that can be approached without conflicting with the PBH data do not depend 
significantly on the shape of primordial perturbation spectrum. The constraints on the GW back- 
ground obtained in previous works are extended to a wider GW frequency range. We discuss the 
applicability of the currently available pulsar timing limits for obtaining the constraints on scalar 
power spectrum and PBH abundance and show that they can be used for strongly constraining the 
PBH number density in the PBH mass range ~ (0.03 — 1O)M0. 
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I. INTRODUCTION 

It is well known now that gravitational waves (GWs) 
can be effectively generated by density perturbations dur- 
ing the radiation dominated era. Tensor and scalar per- 
turbations are decoupled at the first order, but it is not 
so in higher orders of cosmological perturbation theory. 
Namely, the primordial density perturbations and the as- 
sociated scalar metric perturbations generate a cosmo- 
logical background of GWs at second order through a 
coupling of modes 1-3]. In particular, a second order 

(2) 

contribution to the tensor mode, , depends quadrat- 
ically on the first order scalar metric perturbation, i.e., 
the observed scalar spectrum sources the generation of 
secondary tensor modes. By other words, the stochastic 
spectrum of second order GWs is induced by the first or- 
der scalar perturbations. Calculations of flaw at second 
order and discussions on perspectives of measurements of 
the second order GWs are contained in works [^-Q ■ 

It is natural to conjecture that the detection of GWs 
from primordial density perturbations on small scales 
(not directly probed by observations) could be used to 
constrain overdensities on these scales, in a close analogy 
with the case of primordial black holes (PBHs). How- 
ever, at the present time, gravitational wave background 
(GWB) is not yet detected. So, on the contrary, one can 
constrain GWB using existing limits on amplitudes of 
primordial density perturbations. Such limits are avail- 
able, in particular, from studies of primordial black hole 
production in the radiation era. 

It is generally known that PBHs form from the density 
perturbations, induced by quantum vacuum fluctuations 
during inflationary expansion. For an efficient produc- 
tion of PBHs in the early Universe [lol - [l3| the spectrum 
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of the density perturbations set down by inflation must 
be "blue", i.e., it must have more power on small scales. 
This implies that the spectral index of the scalar per- 
turbations must be larger than 1, in strong contradiction 
with the latest WMAP results [l3 - [l6j . Such a conclusion 
is correct, however, only in rather special case: namely, 
it is based on the prediction of slow-roll single field infla- 
tionary scenario, according to which the power spectrum 
of curvature perturbations is nearly scale- invariant, i.e., 
the spectral index n is close to unity and the variation in 
the spectral index dn/dlogk is small. 

Although a prediction of the approximate scale in- 
variance of the primordial power spectrum is a neces- 
sary requirement to any inflationary model, some de- 
viations from pure scale invariance are consistent with 
the observational data. These deviations are described 
by adding localized features to the primordial spectrum 
(see, e.g., [17] and references therein) and/or by intro- 
ducing spectral features modifying a single power law. 
Models with such peculiarities (sometimes called broken- 
scale-invariant (BSI) models) were proposed, in main as- 
pects, in eighties fl8t - i23i] . Such models generally include, 
in addition to the usual inflaton field, other scalar fields 
driving successive stages of inflation and triggering phase 
transitions. 

Evidently, the BSI models of inflation could predict, 
generically, the essential production of primordial black 
holes at small and medium scales. In particular, second 
order phase transitions during inflationary expansion had 
been first considered in j20l . [23 j in models with two scalar 
fields. In scenarios of such type, during a short stage, 
corresponding to the beginning of a phase transition, the 
mass of the trigger field becomes negative, and adiabatic 
perturbations are exponentially amplified resulting in the 
formation of a narrow spike in the primordial spectrum, 
and, as a consequence, in a copious production of PBHs 
(2^ [25! . There are many multiple field scenarios predict- 
ing the existence of spike- or bumplike features in the 
primordial spectrum (e.g., supersymmetric double hybrid 



2 



models [26|, multiple inflation models based on super- 
gravity [23, etc). Some of these models are spec ially 
constructed to predict efficient PBH production [23, [23| . 

An existence of the narrow spikes in the primordial 
spectrum is possible not only in multiple field inflation- 
ary scenarios. Such a feature can, in p rinciple, exist even 
in single field models (see, e.g., ^fllHl)- If, e.g., the 
inflationary potential has an unstable maximum at ori- 
gin (e.g., the double-well potential) then, with some fine- 
tuning of parameters and initial conditions, the inflation 
process may have two stages, with a temporary stay at 
the maximum, that may lead to the corresponding peak 
in the primordial spectrum and, depending on the am- 
plitude of the peak, to the PBH production. 

The details of the PBH formation from the density 
perturbations had been studied in [H, H^l, the astro- 
physical and cosmological constraints on the PBH den- 
sity had been obtained in many subsequent works (see, 
e.g., the recent reviews [I^, Hj]). The order of magni- 
tude of the corresponding constraint on the value of the 
density perturbation amplitude is well known |35], but, if 
the primordial spectrum contains the peak-like feature, 
the concrete value of the PBH constraint clearly depends 
on the parameters characterizing the form of this fea- 
ture (in particular, on the width of the peak). Such an 
information may be rather useful for the model makers. 

The aim of the present work is two-fold. In the first 
part of the work we obtain constraints on a power spec- 
trum of the primordial fluctuations (for the particular 
case, when the spectrum has a peak feature), for a wide 
range of PBH masses (10^ 10'^* g). Recently, the con- 
straints on the curvature perturbation from PBHs had 
been compiled and updated in Ref. |36|]. Authors of [3^ 
assume that the PBHs form at a single epoch and that, 
over the scales probed by a specific PBH abundance con- 
straint, the curvature power spectrum can be written as 
a power law (with a spectral index close to 1). In con- 
trast with this, we assume that the curvature perturba- 
tion spectrum has a peak, and the position of this peak 
determines the epoch of the PBH production. A width 
of this peak is a model parameter, and the peak value 
is constrained by corresponding data (on nucleosynthe- 
sis, photon extragalactic background, cosmological en- 
ergy density parameter). 

In the second part of the work we use the constraints 
on the curvature perturbation derived in such a way 
for constraining the energy density of the induced GW 
background (different values of PBH masses correspond 
to different values of a frequency of this background). 
In our previous work (37j a part of these constraints 
were obtained, for rather narrow range of frequencies 
(~ I0~'^ 10'^ Hz), whereas in this paper we do it for 
the interval - IO"i° ^ 10* Hz. 

The plan of the paper is as follows. In Sec. [Tl] we 
present the main relations which are necessary for a PBH 
mass spectrum calculation. In Sec. IIIII we introduce 
our parametrization of the power spectrum of primor- 
dial curvature perturbations, having a peak feature, and 



demonstrate a dependence of the PBH mass spectrum 
on a width of the peak. In Sec. IIVI we obtain the con- 
straints on the peak value of the primordial curvature 
spectrum from non-observation of PBHs and products of 
their Hawking evaporation. In Sec. |V]we give the main 
formulas used for the calculation of the induced GWB. 
Constraints on ilciv derived from PBH constraints on 
the primordial curvature spectrum are presented in Sec. 
IVIl The last Section contains our conclusions and dis- 
cussions. 



II. PBH MASS SPECTRUM CALCULATION 

The calculation of PBH mass spectrum in Press- 
Schechter formalism [38] is based on the expressions [si" 



nBH{MBH)dMBH = 
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(1) 



(2) 



Here, the following notations are used: 5b. is the initial 
(at the moment ti) density contrast smoothed on the co- 
moving scale i?, M is the smoothing mass (initial mass 
of the fiuctuation corresponding to the scale i?), (Tfl (M) 
is the mean square deviation (the mass variance). 
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Vs{k)W\kR)-, 
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Vsik) is the power spectrum of primordial density per- 
turbations, W(kR) is the Fourier transform of the win- 
dow function (in this work we use the gaussian one, 
W{kR) = exp(-fc2i?2/2) ), is the initial energy den- 
sity. It is assumed that the process of reheating is very 
short in time, so, the end of inflation practically coincides 
with a start (at t — ti) of the radiation era. 

Fourier transform of the (comoving) density contrast 

is 



5k{t) 



2 / fc 



*fc(t), 



(4) 



where ^'fe(i) is the Fourier transform of the Bardeen po- 
tential. Here, we explicitly take into account the time 
dependence of the Bardeen potential. 

The power spectrum of the density perturbations, cal- 
culated at some moment of time, is 



Vsik,t) 



1 2 



r^ik.t), 



(5) 



where r is the conformal time (t = {aH) ^ for the radi- 
ation epoch). 
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The comoving smoothing scale, R = l/ka, is con- 
nected with the smoothing mass M by the expression 



/ M \ -2/3 



{aiHi 



(6) 



where M^, Ui and Hi are the horizon mass, cosmic scale 
factor and Hubble parameter at the moment ti. 

In the approximation of instantaneous transition from 
inflationary era to the radiation dominated epoch, the 
connection between density perturbation at any time and 
curvature perturbation at initial moment of time U is j42j | 



Vs{k,t) 



n 2 



where the expression for ^fe(T) is given by [42 
27^fc(r,) , 



*fe(T) 



-[{x — Xi) cos(2; — Xi) — 

-(1 + xxi) sm{x — Xi)]. 
kr kTi 



(7) 



(8) 



v/3' 



Here, TZk{Ti) is the Fourier component of the curvature 
perturbation on the comoving hypersurfaces, at the end 
of inflation (see, e.g., jSsj). 

The connection between values of the smoothing mass 
M, density contrast and PBH mass Mbh can be ex- 
pressed in the general form 



(9) 



The concrete expression for the function ip depends on 
the model of the gravitational collapse. In the model of 
the standard spherically-symmetric collapse the connec- 
tion is quite simple: 



BH 



(10) 



Here, Mh is the horizon mass at the moment of time, t = 
t/i, when regions of the comoving size R and smoothing 
mass M cross horizon. According to Carr and Hawking 
1/3 = 5th < -^f < 1. The derivation of Eq. ^ 
is given in the Appendix of [l^. From pil)) . using the 

relation Mh = M^^^ M'^/^ , one has the expression for the 
function for the Carr-Hawking collapse: 
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FIG. 1: The calculation of PBH mass spectra for a peaked 
■P7j(fc)-spectrum, for = 10'^'^ g and different values of E. 
The value of was chosen so that for all cases Q,pbh ~ 
0.024 (for E = 1, -P^ = 0.0488; for E = 3, = 0.0294; For 
E = 5, -P^ = 0.0267; For E = 7, = 0.0258.) The dashed 
curve shows, for comparison, the form of the classical Carr's 
Af^^''^-mass spectrum 32]. 



maximum of PBH mass spectrum is still around the hori- 
zon mass corresponding to the maximum in primordial 
power spectrum, but the PBH mass spectrum also has a 
"tail" of small masses. In this work we will use classical 
collapse model (Eq. ([TT|) ). However, in the end of the 
paper we will explore the dependence of the results on 
5th - the density contrast threshold of PBH formation, 
which gives the main uncertainty due to the exponential 
sensitivity of PBH abundance to it. 



III. PRIMORDIAL POWER SPECTRUM WITH 
MAXIMUM 

It is convenient to use some kind of parametrization 
to model the realistic peaked power spectrum of finite 
width. We use the distribution of the form 



</p(M, 5^- Mi) = [5^)^IHI^/^mI'\ (11) lg^7j(fc) = B + {\gVl- B) cxp 



{^k/k^r^ 



(13) 



In the picture of the critical collapse |4£ 
responding function is 



14611 the cor- 



(^(M,<5f ;M,) = fc,(5f - 5,y^M''I^M\l\ (12) 

where 5c^ "lc and kc are model parameters. The mass 
spectrum of PBHs for the critical collapse model has 
been calculated, e.g., in It was shown that for the 
primordial scalar perturbation spectrum with a peak the 



Here, B w —8.6, V\ characterizes the height of the peak, 
fco is the position of the maximum and S is the peak's 
width. Parameters of such a distribution have been con- 
strained in the previous work of authors [42J from non- 
observation of PBHs and products of their Hawking evap- 
oration (photons and neutrinos). 

In many cases limits on PBH abundance have been 
obtained using the assumption that PBH mass distri- 
bution is close to the (5- function form. This does not. 
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however, mean that the original spectrum Vn{k) should 
be very narrow. For example, in the calculation of PBH 
mass spectra given in Fig. [H all distributions have a 
distinct maximum (even for rather wide peak, S = 7), 
which shifts to smaller Mbh with the growth of E. In 
this sense, if our goal was just to explore the PBH abun- 
dance, we could assume that T'Tz{k) ^ (5(fc — fco) from the 
start, without introducing big mistakes (note, however, 
that with such an approach the shift of the distribution 
maximum would not be noticed). This is, however, not 
an adequate approximation for our study, because we are 
also interested in calculation of the 2-nd order GWB pro- 
duced by the same power spectrum and studying of its 
dependence on the width of the peak. 

One should note that in Fig. [T]we show on the vertical 
axis the quantity ubh x (a^/ao)'^, which is the physical 
(rather than comoving) number density of PBHs, and is 
independent on the reheating temperature Tjih in the 
limit h = a^H^ > fco = k{M^) [42]. 

The fraction of an energy density of the Universe con- 
tained in PBHs today, fipBH, assuming that a mass of 
the produced black hole does not change in time, is 

^PBH = — ( — \ I MBHnBH{MBH)dMBH (14) 

Pc \aoJ J 

{pc is the critical density). This formula is rather accu- 
rate for black holes with initial mass Mbh ^ Af*, where 
w (3toao)^^^ ~ 5 X 10^"* g is the initial mass of PBH 
which reaches its final state of evaporation today (47,], 
ao = 8.42 X 10^^ g'^s"^, and to is the age of the Universe. 



IV. CONSTRAINTS ON 

For constraining we use the existing limits on PBH 
abundance. In the region of Mbh which is of interest 
for us (10^ ^ Mbh ^ 10^* g) these limits can be di- 
vided in three groups: i) constraints on PBHs from big 
bang nucleosynthesis (due to hadron injections by PBHs 
(isj . photodissociation of deuterium [49^] and light nu- 
clei, fragmentations of quarks and gluons evaporated by 
PBHs [Sfl]), 10^ < Mbh ^ 10^^ g, and from infiuence of 
PBH evaporations on the CMB anisotropy, 2.5 x lO^'^ ;$ 
Mbh ^ 2.5 x 10^'' g [11], ii) constraints on PBHs from 
extragalactic photon background, 10^^ < Mbh ^ 10^^ g, 
in) constraints on non-evaporating PBHs (gravitational 
and lensing constraints) . Constraints on PBHs from data 
on extragalactic neutrino background [4l|, |4^ [5l| , in the 
region 10^"'^ ^ Mbh ^ 10^^ g, are somewhat weaker than 
nucleosynthesis constraints. 

We do not consider in this paper the PBH constraints 
from PBH masses smaller than 10^ g because the cor- 
responding GWB frequencies are too high (> 10'' Hz). 
Also, we did not consider in this paper the constraints 
from galactic gamma rays (they could be essential in the 
narrow mass region near Mbh 10^^ g, but strongly 
depend on the unknown clustering factor [SJl]) and po- 



tential constraints [3J] from future measurements of 21 
cm line. 

For a derivation of the constraint on in the re- 
gion 10^ <; Mbh ^ 10^'^ g we use the latest update of 
the nucleosynthesis constraints given in the review [si] ]. 
In this region of PBH masses we, following the tradi- 
tional practice, approximate the initial PBH mass spec- 
trum by (5- function (i.e., we assume that all PBHs have 
the same mass Mbh)- The mass of the PBH formed is 
approximately equal to the horizon mass at horizon en- 
try, Mbh ~ Mh- Correspondingly, the constraints are 
expressed in terms of the function (3{Mbh) which is the 
mass fraction of the energy density of the Universe going 
to PBHs. It is given by the relation 



f3{MBH) - 2 J P{5uor{R))d5hor{R) 

= ertc —= 

\V2aR{MBH) 

where the mass variance is given by [s^ 



oI{Mbh) = y / [kRfjl 



(15) 



(16) 



xexp(-fc^i^^)7'7^(fc) 



dk 



R is the smoothing scale, i.e., the horizon size. 

The order of magnitude of the constraint on /3 is, ac- 
cording to HI: /? < 10-18 for lO^g < Mbh < lO^^g and 
/3 < 10-23 for lO^^g < Mbh < lO^^g. In the narrow 
region near ~ 3 x lO^"^ g there is the strong constraint 
following from CMB anisotropy damping, /3 ^ 10-^*. 

For a derivation of the constraint on Vji in the region 
10" < Mbh ^ 10^'^ g we use the extended PBH mass 
spectrum given in Sec. Hi] and the technique developed by 
authors in 42] (in particular, the condition that photons 
(primary as well as secondary ones [s^ ) evaporated from 
PBHs do not exceed the observed extragalactic gamma 
ray background was used) . Our method is essentially the 
same as used in the pioneering work [ssj . Extragalactic 
photon background data which we used had been ob- 
tained in [54] (EGRET coUaboration) and in [55] (Fermi 
LAT collaboration). Comparison shows that in the en- 
ergy region ^ 0.1 1 GeV (which is of interest for us) 
data of both works are consistent with each other. So, 
for concrete calculations we used the data from [HI] . 

At last, the gravitational constraint is just a condition 
that ^pBH does not exceed the energy density of non- 
baryonic dark matter (in the limiting case, PBHs account 
for all dark matter): 



PBH 



< n 



COM 



0.25. 



(17) 



The value of is obtained from and (fTT]) . 

Lensing constraints arise due to the fact that mi- 
crolensing observations of stars in the Magellanic Clouds 
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FIG. 2: The limits on the value of Vn from PBH non- 
observation, for three different values of E. The constraints 
based on the /3-function approach are shown by thin lines. 



probe the fraction of the Galactic halo in massive com- 
pact halo objects (MACHOs) of sub-solar masses [s^. 
The fraction of the halo in PBHs is 



BH) — 



PBH 



COM 



(18) 



From the analysis of MACHO 57] and EROS ^ mi- 
crolensing surveys, /{Mbh) < 0.1 for 10~^Mq < 
Mbh < Mq and UMbh) < 0.04 for IQ-^Mq < 
Mbh < O.1M0. There are some additional constraints 
on /{Mbh) for other mass ranges, reviewed, e.g., in [s^ . 
but they are generally rather weak {/{Mbh) ~ 1) and 
we do not use them here. 

The resulting constraints are shown in Fig. [2] The con- 
straints based on the I3{Mbh) function are designated 
by thin lines (the corresponding constraints on j3 have 
been taken from Fig. 6 of [3J]). The constraints ob- 
tained by using our technique (which is briefly explained 
is Sec. ini and EGRET data, and gravitational con- 
straints are shown by thick lines. For calculations of lim- 
its on flaw in the following Sections we use, for smaller 
masses (before the crossing points of thin and thick lines) , 
the "/3-constraint" , and for larger masses, after the cross- 
ing points, the constraints given by the thick lines. The 
value of the crossing point is about 4 x 10^'^ g, almost 
independently on the value of E. 

One can see from Fig. [5] that the results for con- 
straint depend on the method of the calculation, in the 
region ~ 10^"* 10^^ g. Namely, the constraints based on 
the /3 function are somewhat weaker. We did not study, 
in the present work, the dependence of this result on the 
form of the PBH mass spectrum. 



The ■p.^-constraints shown in Fig. [2] (for the cases of 
S = 3,5) are close to the ones derived in [sl]. Authors 
of [IBl have assumed that the power spectrum is scale- 
invariant over the (relatively small) range of scales which 
contribute to a given constraint (this is analogous to the 
assumption that primordial spectrum has a rather wide 
peak) . They argued that deviations from scale-invariance 
(which are consistent with a slow-roll inflation hypothe- 
sis) lead to small changes in the constraints. However, in 
the case of a narrow peak (our case of S = 1), which can 
be produced, e.g., by a violent slow- roll violation during 
single-fleld inflation, the limit on Vji will, of course, be 
much weaker, as we see from Fig. [2] (because the total 
power contained in the spectrum depends on its width). 



V. INDUCED GW BACKGROUND 

A. Basic formulas 

According to f6| , the power spectrum of induced GWs 
is given by the expression 

CO 1 

Vh{k,T)^ j dk j d^iV^{\k-k\)V^{k)J'{k,k,fi,T), 

-1 



(19) 



where 

T{k, fc, ^, r) = ^ — r— ^ ^ X 

a^{T) |k-kp 

T 

X / dfi a(fi)gfe(T,fi)/(k, k,fi) X 



(20) 



X J df2a{f2)gk{T,T2) /(k,k, fs) -l-/(k,k-k, fs) 

To 

and 

/(k,k,r) = 12vI/(fcT)*(|k-k|r)+ (21) 
+&T-9{kT)^'{\k - k|r) + 4:T^^'{kT)-9'{\k - k|r). 

In Eqs. [201 1^1]) the following notations are used. 

V<Sf{k) is the power spectrum of the Bardeen potential, 
defined at some moment of time r = r^' near the begin- 
ning of the RD stage (by definition, it is the primordial 
spectrum) , 



2^2 

k 



3 5^{k + k')V^{k), 



is the Fourier component of 5* 
1 



*(x) = 



(27r)3/2 



(22) 



(23) 



/i = k-k/(fcfc) is the cosine of the angle between the 
vectors k and k. The power spectrum of GWs is defined 
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by the standard expression 



1 

2^ 



(24) 



where h]s_{T) is the Fourier component of the tensor met- 
ric perturbation, 



htj{x, t) = 



(27r)3/2 



(25) 

eij (k) and (k) are two polarization tensors correspond- 
ing to the wave number k. 

The function / in Eq. contains transfer functions 
^'(fcr), which are defined by 



*(/cr) 



(26) 



where = ^'^(t/) is the initial (primordial) value of the 
potential. During RD epoch, the solution for the Bardeen 
potential, having the initial condition 4'fe(Ti) = 0, where 
Ti is the moment of the end of inflation which is close to t[ 
(but Ti < t[), is given by Eq. ([SJ. The value of the poten- 
tial at Ti is chosen to be zero because ^fe is typically very 
small during inflation (isj and it is a continuous function 
during the transition from inflationary to RD stage (we 
assume, for simplicity, that the reheating is instant). We 
have chosen, for the numerical calculation, the value of 
r/ using the condition \g{T-/Ti) = 0.05 and thereby ne- 
glect the formation of PBHs and induced GWs in the 
interval of time from Ti to r^'. In this work we are inter- 
ested only in wave numbers k <C kend, for which there is 
no dependence of the results on Ti (because the pertur- 
bation amplitudes, such as 4'fc(T), have enough time to 
reach their asymptotic limit before horizon re-entry for 
each mode). For an example of the case where values of 
k ~ kend are important, see, e.g., Ref. [1^ (where the 
case of the running mass model is considered). 

The function gk{T,f) in Eq. (f20l) is the Green func- 
tion which depends on the cosmological epoch. For RD 
Universe, 



to calculate flaw for the moment of time Tcaic ^ k~^ at 
which the mode is well inside the horizon, and is freely 
propagating. We can then easily relate energy densities 
of GWs at different times with simple calculation, us- 



ing the fact that hk 



far inside the horizon (this 



approach works well only for rather large values of fc. 



k > kr 



lOO/Ce 



IMpc , but we are only interested 



in such large wave numbers here). 

The final expression for QQ^r{k) is 

n%^ik) = 2nn ( X i^!^7',(fc,reaic 



9* 



12 



(29) 

This formula gives the correct energy density, accurate 
to the oscillations in it. In practice, Tcaic can be either 
fixed or dependent on fc, e.g., for the last case. 



sub 



100. 



(30) 



It proves to be more convenient to use the "randomized" 
value of iVsub, i-e.. 



"^calc 



(Nsuh + N,nd) ■ k- 



(31) 



where Nsnh is constant and iYmd is a random number in 
the interval [0, 27r] calculated independently for every fc. 
In this case the result of the calculation is a stochastically 
oscillating function whose envelope always can be easily 
found, and it is the envelope that we are interested in. 
The exact shape of the function will, actually, depend on 
the choice of Tcaic ^ k~^, and the larger Tcaic we take, the 
more frequent are the oscillations, but the envelope which 
we are interested in does not change. This was explicitly 
shown in work Q (see, in particular. Figure 2 from it, 
which shows the same GW spectrum calculated using 
approach of Eq. pip and the one using Tcaic = const. 
It is seen from that figure that the resulting spectrum is 
the same). 



B. Connection between frequency and horizon 
mass 



k{T,T) = ^sin[fc(T- f)] 

K 



T < Teq. (27) 



The energy density of GWs per logarithmic interval of 
k in units of the critical density is given by 



1 



^Gwik,T) = — 



12 \a{T)H{T) 



Vh{k,T). (28) 



Here, the power spectrum of GWs, VhikjT), is obtained 
from the formula ()19|) . However, for very large wave num- 
bers fc which we are interested in, the direct use of (fT9|) 
will require numerical integration for functions having a 
huge number of oscillations (e.g., for fc ^ lO^^Mpc"^ this 
is about ~ fcTo ^ 10^" oscillations). This is hard to do 
numerically. Fortunately, we do not have to do integra- 
tion until the present day. As discussed in fol] , it is enough 



For a wave with comoving wave number k and wave- 
length A — 27r/fc, propagating at the speed of light c, the 
corresponding frequency is / = c/A, or 

ck / fc \ 

/ = — 1.54 X 10"^^ ^ Hz. (32) 

27r VMpc^V 

From the constancy of the entropy in the comoving vol- 
ume, we have the relation between the scale factor a, 
temperature T and the effective number of degrees of 
freedom : 



a ^ g^, ' 1 

From the Friedmann equation {H^ ^ p), we have 

Tj -2 -1/6 

H ^ a g* , 



(33) 



(34) 
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and the horizon mass corresponding to the scale factor a 
evolves during the RD epoch as 



(35) 



From and ([351) . the wave number of the mode en- 
tering horizon at the moment of time t (at this time, 
k = aH) is related to the horizon mass at the same mo- 
ment of time by 



k — k^ 



■eg 



, , X -1/2 , 

Mh \ I 9* 



9*eq 



M, 



eq 



-1/12 



2x 1023(M;,[g])-i/2 Mpc-\ 



(36) 



where in the last equality we have adopted that g^eg ~ 3, 
w 100, 



Meq = 1.3 X 10" 



50„ 



-1/2 



X 10 



0.0095 Mpc" 



(37) 



(38) 



The factor is a function of cosmic temperature T 
and in this work we assume that there are no degrees 
of freedom beyond the Standard Model, so g^, « 100 for 
T > 100 MeV, « 10 for 100 MeV > T > 1 MeV and 
5* « g^:f,g « 3 for T < 1 MeV. The connection between 
horizon mass and T is 



Mh « 3 X IO^'Mq ( — 



-1/2 



T 



1 MeV 



(39) 



and we can estimate that M/i(l MeV) ~ 10^ Mq, and 
M,,(100 MeV) « 5Mq. 

The frequency of the wave corresponding to the wave 
number k can be related to the horizon mass by the re- 
lation following from and 



/«3x 10« Hzx (M;,[g]) 



^1/2. 



9 X lO^^g 



(40) 



For scalar-induced GWs, the single mode in scalar spec- 
trum does not correspond to the only one mode in Vh- 
For example, for the (5-function-like spectrum Viz{k) ^ 
d{k — ko), the GW spectrum is continuous and stretches 
from to 2fco [51] . However, the order of magnitude of 
wave numbers of induced GWs is the same as of scalar 
perturbations, so (|40p gives an estimate of GW frequency 
that will be generated from perturbations entering hori- 
zon at its mass scale Mh- Furthermore, if PBHs form 
from a scalar spectrum of perturbations at a horizon mass 
scale Mh, the typical PBH mass will be of order of Alh, so 
relates the typical PBH mass with the characteristic 
frequency of second-order GWs produced. 



logio f [Hz] 



o 



o 
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VI. CONSTRAINTS ON ilr 



FROM PBHS 



13.5 14.0 14.5 15.0 15.5 16.0 
log 10 k [Mpc"^] 

FIG. 3: a) The calculation of the induced GWB correspond- 
ing to primordial power spectra and PBH mass spectra given 
in Fig. [1] (cases of E = 1,3,5). The pulsar timing limit for 
Qgw is also shown, b) The calculation of the induced GWB 
for = lO^'^ g (for all three cases, ilpBH 0.24, the cor- 
responding sets of parameters are = 0.028 for S = 1, 
= 0.0172 for E = 3, and = 0.0159 for E = 5). 



that for a narrow peak, the distribution looks much like 
the doubl&jjeaked one produced by a (5-function power 
spectrum [5|, l9| . The shape is smoothing with the growth 
of S (and it will be a scale-invariant spectrum for the 
scale-invariant input The value of Hqw in the case 
of rather wide peak (E ^ 1) is proportional to [Vfif' 
and can be estimated as 



In Fig. [3] we show the result of View calculation for 
the finite- width distribution of the form ([T5)) . It is seen 



r!Gvy(fc>fcc,ro) = 0.002 

9* 



1/3 



■V- 



(41) 
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It is seen from Fig. [3] that the maximal values of 2- 
nd order GWB amplitudes that can be reached do not 
depend significantly on the width of the primordial power 
spectrum. The growth of flaw with increase of E, which 
is naturally expected (see [9,]), is partly compensated by 
the simultaneous decrease of the curvature perturbation 
amplitude (really, Fig. 2 shows that constraint on 
Vtz decreases with the growth of S). This allows to put 
a constraint on the value of ^gW: which is independent 
on the width. Such constraint is shown in Fig. 2] for the 
whole range of GW frequencies considered, ~ 10~^° -;- 
10^ Hz. 

The shaded triangle in Figs. |3^,|4]designates the pulsar 
timing limit obtained in [59|, 

^owif)h^<^.S.10-^[^;^^J^)\ (42) 

for / > 4.4 X 10-^ Hz at 90% C.L. It is seen that PBH 
mass spectra from Fig. [I](M° = 0.5Mq) are inconsistent 
with this limit: the same scalar perturbation spectrum 
overproduces GWs. We have estimated that the limit 
(|42|) excludes the significant amount of PBHs in the re- 
gion of masses ~ (0.03-^ lO)Af0 due to overproduction of 
2-nd order GWs. Note that this conclusion differs from 
the earlier result of 7] who associate the pulsar timing 
limit (|42[) with PBHs of mass ^ lO^M©. This is obviously 
caused by the definition / = 2TTck used by these authors, 
instead of the usual one, given by Eq. ((32|) . From our re- 
sults it follows that the primordial origin of intermediate 
mass black holes (IMBHs), with masses {10'^~10'^)Mq, 
is not excluded by pulsar timing limits. 

PBHs with mass close to « 10^^ g can be responsible 
for cosmic dark matter, and, as shown in [i^,!!^], if they 
are clustered in the Galactic Center, they can explain 
the 511-keV photon line observed from its direction (such 
photons, in this case, are produced by the annihilation of 
positrons evaporated from PBHs [6l|,[62l)- It is seen, in 
particular, from Fig. [3Jd that to explain this phenomenon 
with clustered PBHs, flew must be approaching a value 
of « 10~^ near / ~ 1 Hz. This region will be probed in 
future by BBO experiment. 

Currently, the analysis of pulsar timing data is the 
only experiment which allows to set a stronger bound 
on primordial spectrum than PBHs (see also [63|). For 
comparison we have also shown the limit obtained by 
the ground-based interferometer LIGO during its fifth 
science run (S5) [63 |. 

Qgw < 6.9 X 10"^ (43) 

This limit applies to a scale-invariant GW spectrum in 
the fre quen cy range 41.5 — 169.25 Hz. The target sen- 
sitivity [65[ of the planned Advanced LIGO experiment, 
^Gw ^ 10^^ 10~^, is also shown in Fig. |31 It is seen 
that in future the obtained limit for GWB will be experi- 
mentally approached in this experiment, and also in most 
other regions of the frequency range considered, by other 
experiments such as LISA [66j], Big Bang Observer (BBO, 



logio f [Hz] 



-9. -7. -5. -3. -1. 1. 3. 




4 6 8 10 12 14 16 18 



logiok[Mpc 1] 

FIG. 4: The limits on the 2-nd order GWB from primor- 
dial black holes, obtained in this paper. Solid thin line shows 
the result assuming PBH formation threshold is 5th = 1/3. 
Dashed line corresponds to the assumption 5th ~ 0.45. Thin 
line shows the maximum values of GWB that can be reached 
for ^pBH ~ 10""" (and 5th = 1/3). Also shown are current 
pulsar timing limit, LIGO S5 limit, Advanced LIGO planned 
sensitivity to Qgw, and frequency ranges in which other fu- 
ture experiments (LISA, BBO, SKA) will operate. The esti- 
mate sensitivity of all future experiments, in the correspond- 
ing frequency ranges, is much better than the PBH bound 
shown. 

see, e.g., [S^l), Square-Kilometer- Array (SKA, see, e.g., 

[si). 

To show the uncertainty in PBH constraints, we have 
also plotted in Fig. Uthe resulting constraints obtained 
assuming a somewhat larger PBH formation thresh- 
old, Sth = 0.45, and maximal GWB corresponding to 
^PBH = 10-5 (and Sth = 1/3). 



VII. SUMMARY AND DISCUSSION 

We have performed calculations assuming that PBHs 
and 2-nd order GWs are produced at times not very close 
to the beginning of the radiation era U. If production 
of GWs and PBHs takes place near ti (i.e., right after 
the end of infiation) , the limits derived in our paper can 
potentially be altered. One such example is a particu- 
lar case of the running mass inflation model (proposed 
in [6^, [tO] and further studied in many papers (see [3l| 
and references therein)), which predicts a rather strong 
scale dependence of the spectral index, possibly leading 
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to largest values of Vu{k) just near ki (« kend)- The 
analysis performed in [2], however, shows that the max- 
imum values of GWB amplitude that can be reached in 
this case are also very close to maximum values derived 
in this paper. 

In summary, we have performed simultaneous calcula- 
tions of PBH mass spectra and induced GW background, 
obtaining the constraints on values of VTz{k) and D,Gw{k) 
from known limits on the PBH concentration in various 
cosmological scale ranges. We have explored the depen- 
dence of these limits on the shape of the primordial spec- 
trum (in particular, on its width). It was shown that 
though constraints on the peak width may significantly 
depend on the shape of the spectrum, the maximal pos- 
sible values of GWB are almost unsensitive to it. This 
allowed us to place quite model-independent limits on 



induced GWB in the wide frequency range. We have 
discussed the applicability of currently available exper- 
imental data, in particular, pulsar timing limits, to the 
constraining of PBH abundance. We have shown that 
the primordial origin of IMBHs is not forbidden by the 
pulsar timing limits. 

Comparing our results with the previous results of [36j 
and 0, one can see that our constraints for 7-"^ almost 
coincide with those of Ref. [s^ in cases when the width 
of the peak in primordial spectrum is large, E = 3, 5 
(see Fig. [5]); our constraints for il.GW are systematically 
weaker than those of Ref. |8i] . Also, in our case the depen- 
dence of ilcvi^-constraint on the shape of the primordial 
spectrum is much smaller than the analogous dependence 
in Ref. 
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